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With the ever increases in the capabilities of computers for numerical computations, we 
are on the verge of using these tools to model manufacturing processes for improving 
the efficiency of these processes as well as the quality of the products. One such pro- 
cess is casting for the production of metals. However, in order to model metal casting 
processes in a meaningful way it is essential to have the basic properties of these mate- 
rials in their molten state, solid state as well as in the mixed state of solid and liquid. 
Some of the properties needed may be considered as intrinsic such as the density, heat 
capacity or enthalpy of freezing of a pure metal, while others are not. For instance, the 
enthalpy of solidification of an alloy is not a defined thermodynamic quantity. It’s value 
depends on the micro-segregation of the phases during the course of solidification. The 
objective of the present study is to present a thermodynamic approach to obtain some of 
the intrinsic properties and combining thermodynamics with kinetic models to estimate 
such quantities as the enthalpy of solidification of an alloy. 

A schematic diagram illustrating our general approach in obtaining the thermodynamic 
description of an alloy is shown in Figure 1 . The starting point is to have suitable and 
appropriate thermodynamic models for all the phases involved. All of these models have 
unknown parameters. The parameter values must be obtained from experimental ther- 
modynamic data. Although most thermodynamic models are empirical or semi-empirical, 
we attempt to use models with some physical basis and certainly with the least number 
of parameters to describe such a phase. Once we have the thermodynamic properties of 
all the phases in a binary described, we can then calculate stable binary phase equilibria 
and compare them with the measured diagrams. If the thermodynamic data and phase 
diagram data are correct and the models used are suitable, the calculated diagram 
should agree with the measured equilibria. Since all experimental measurements have 
uncertainties, it is often necessary to use thermodynamic and phase equilibrium data 
simultaneously to optimize the parameter values of the thermodynamic models. Of 
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Figure 1 

course, we always use the thermodynamic values as the starting point for our optimiza- 
tion. We can also calculate phase equilibria at either very high temperatures where mea- 
surements are difficult to carry out experimentally or at very low temperatures where the 
kinetics of reaction are so slow that achievement of equilibrium is difficult if not impossi- 
ble. We can also calculate metastable equilibria from the thermodynamic models particu- 
larly those involving a supercooled liquid phase. A knowledge of the metastable equilibria 
is important in controlling the resulting microstructure of alloys during solidification. In sev- 
eral systems, our calculated metastable equilibria are in accord with experimental mea- 
surements. In addition, on the basis of the binaries, we may extend the calculations/ 
predictions to ternary and higher order systems. This is particularly important since most 
commercial alloys are of multi-components. 

This type of approach allows us to have the thermodynamic descriptions of the phases 
involved as well as to obtain the values for the distribution coefficients as a function of 
composition and temperature. The distribution coefficient is a basic parameter in model- 
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ing the microsegregation of an alloy. Yet it is usually not available from experimental 
measurement for a multi-component alloy as a function of temperature. Coupling the 
thermodynamic descriptions of the phases with a kinetic model, we can compute such 
important quantities such as the enthalpy of solidification which is needed for macro- 
modeling of solidification. Moreover, coupling of the thermodynamic descriptions of the 
phases with semi-theoretical models, we can estimate/calculate the surface tensions of 
molten alloys provided the values of the component elements are known. 
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Nothing is better 
than 

good experimental data 


But We Have To Be 
Realistic 


We don't have either the time or 
the resource to measure all the 
needed materials properties. 


We must develop 
methodologies 
to estimate these 
materials properties 
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In fact, a combined thermodynamic 
modeling/experimental approach is 
the only sensible way to build up a 
material data base for engineering 
applications! 


Solid-liquid equilibria 

(distribution coefficient) 
Enthalpy of fusion 
Enthalpy of Solidification 
Heat capacity 
Surface tension 


Thermodynamic Approach 
to Obtain Phase 
Equilibrium Data 



Thermodynamic 

models: 

Parameters 


Experimental 
measurements of 
Gibbs energies 


Calculation of 
stable equilibria 


Experimental 
determinations of 
phase boundaries 


Prediction/ 
calculation of 
higher order 
phase equilibria 


Prediction of 
metastable equilibria 
at low temperature 
and stable equilibria 
at very high temperature 


Distribution Coefficient 
( 4 path of solidification) 
Enthalpy of Fusion 
Enthalpy of Solidification 
Heat Capacity 
Surface Tension 
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Thermodyna: 

Models 


U 







Thermodynamic Model (Equation) for Solutions (both Solid and 
Liquid) with Small Deviations from Regular Solution Behavior. 

XS (K 2. 

A G = x(l-x) [u + vx + wx j 

with u, v, w = A + BT + CT InT, where A, B, C are constants. 
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Intermidiate Phases: Variation of the 
Wa^-Schottky 
(Bragg-WStiarns ) Model. 


The Liquid Phase : Associated Solution 
Modd, 

Sublattice Model 
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An Example 


Al-Cu 

Al-Mg 

Mg-Cu 

Al-Mg-Cu 


Al-Cu 

Kaufman and Nesor, 1978 
Murray, 1985 
Chuang, 1984 

Chen, Chuang, Chang and Chu, 1991 



185 












Enthalpy of Solidification 


Case I: Equilibrium Solidification 



B 


Before Solidification: 


AH l (T l , c 0 ) 


After Solidification: 



AH s (T s , C 0 ) 




Enthalpy of Solidification: 

AH = AH S (T S , C 0 ) - AH l (T l , C 0 ) 
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Case U: Eutectic Solidification 




Before Solidification: 



AH l (T e , C e ) 


After Solidification: 



AH a (T E , C^) 

ahP(t e , Cp m ) 


Enthalpy of Solidification: 

AH = [a AH«(T e , C am ) + b AHP(T E , Cp m )] 
- AH l (T e , C e )] 
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Alloy 1 ( 1 . 35Mg , 3 . 98Cu wt%, 1.41 C/sec) 



Enthalpy of 
Solidification 
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Experimental Liquidus Projection of 
Al-Mg-Cu 



Calculated Isotherm of 
Al-Mg-Cu at 400 °C 
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Calculated Liquidus Projection for the 
Al-rich comer of Al-Mg-Cu 


Cu 



Demonstrate the Importance 
of Using Correct 
Distribution Coefficient, k 
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Comparison between the model proposed in 
this study and that by Roozs and Exner 



Model of Roozs 
& Exner 

Present model 

Mass balance: 

Schell equation with 
solid state diffusion 

yes 

yes 

Distribution 

coefficient 

Taken from binaries 
directly as constant 

Computed from 
thermodynamic model 

Heat balance 

yes 

Assumed constant 
cooling rate 

Dendritic 

coarsening 

Described by semi- 
emplrlcal equation 

no 


Alloy 1 ( 1 . 35Mg , 3 . 98Cu wt%) 
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Caselll: Scheil Equation 



kC o 

Scheil Equation: 

C s * = kC 0 (l-^)Ck-l) 



Before Solidification: After Solidification: 



Enthalpy of Solidification: 

AH = (AH° f a + AH E fe) - AH L (T u 
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Casein: Scheil Equation 



A^oCo^ Ce c pmB A^ C ° Cam Q c pm B 

'kC 0 V kC 0 


Scheil Equation: 

C S * = kC 0 (l-fg)(k"l) 


Before Solidification: After Solidification: 



Enthalpy of Solidification: 

AH = (AH a f« + AH® f E ) - AH L (T D CJ 



Caselll: Scheil Equation 



Scheil Equation: 

C s * = kC 0 (l-Js)(k-l) 



Before Solidification: After Solidification: 



AH L (T L > Co) AH E = a AH a (T B C^ + b AH P (T B C pm ) 

a r*i a 

AH =[ AH“(Te. CJ df a ] / f a 
/ 0 

Enthalpy of Solidification: 


AH = (AH® f“ + AH E f E )-AH L (T L , C 0 ) 
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Estimation of 

Cp 


Yeum, Speiser & Poirier (1989) 

Data Needed to Estimate Surface 

Tension (o) of Binary Liquid Alloys 

OA'- Surface tension of the pure 
component A 

03: Surface tension of the pure 

component B 

PA' Density of the pure 
component A 

pg: Density of the pure 

component B 

A^G: Excess Gibbs energy of the 
liquid alloy 


Surface 

Tension 
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10 Binary Liquid Alloy Systems 
were Evaluated 


+ deviations from ideal behavior 
- deviations from ideal behavior 
+,- deviations from ideal behavior 



Surface tensions of Bi-Ag alloys at 1233 K; o « Metzger, 
ideal solution; and *■ calculated. 
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Surface tensions of Ni-Si alloys at 1823 K; o « Popel et al , 
ideal solution; and * calculated. 


Surface tensions of Fe-Cu alloys at 1823 K. o = Dyson; 
ideal solution; and ^calculated. 










Conclusions 


1. A careful assessment of the 
thermodynamic and phase 
equilibrium data of binary and 
higher order systems is essential 
for developing a data base for 
engineering applications. 

Distribution Coefficients 
Paths of Solidification 
Enthalpy of Fusion 
Enthalpy of Solidification 
Heat Capacity 
Surface Tension 

2 A combined modelling / 

experimental approach is the best 
(perhaps the only realistic) 
approach to obtain a 
thermodynamic description of 
multi-component alloy systems. 
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